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Abstract. We study the emission of neutrinos, resulting 
from the scattering of electrons off magnetic flux tubes 
(fluxoids) in the neutron star cores with superfluid (su- 
perconducting) protons. In the absence of proton super- 
fluidity (T > Tcp), this process transforms into the well 
known electron synchrotron emission of neutrino pairs in 
a locally uniform magnetic field B, with the neutrino 
energy loss rate Q proportional to B^T^ . For temper- 
atures T not much below Tcp, the synchrotron regime 
{Q oc T^) persists and the emissivity Q can be ampli- 
fied by several orders of magnitude due to the appearance 
of the fluxoids and associated enhancement of the field 
within them. For lower T, the synchrotron regime trans- 
forms into the bremsstrahlung regime {Q oc T^) similar 
to the ordinary neutrino-pair bremsstrahlung of electrons 
which scatter off atomic nuclei. We calculate Q numer- 
ically and represent our results through a suitable ana- 
lytic fit. In addition, we estimate the emissivities of two 
other neutrino-production mechanisms which are usually 
neglected - neutrino-pair bremsstrahlung processes due to 
electron-proton and electron-electron collisions. We show 
that the electron-fluxoid and electron-electron scattering 
can provide the main neutrino production mechanisms in 
the neutron star cores with highly superfluid protons and 
neutrons at T < 5 x 10^ K. The electron-fluxoid scattering 
is significant if the initial, locally uniform magnetic field 
B > 1013 G. 



1. Introduction 

Neutrino energy losses play a decisive role in cooling of 
young neutron stars (NSs). During first 10^ — 10^ yrs af- 
ter the NS birth, the cooling is dominated by the neutrino 
emission from the NS cores, of density exceeding the nor- 
mal nuclear matter density, po — 2.7 x 10^'^ g cm~'^ (corre- 
sponding to the nucleon number density uq = 0.16 fm^"^). 
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The emission of neutrinos results from the weak in- 
teraction processes involving baryons and/or leptons in 
the NS cores. The presence of the strong internal mag- 
netic field allows an additional process — the neutrino 
synchrotron radiation of electrons. This process was con- 
sidered by many authors (e.g, Kaminker et al. 1991, 1992, 
and references therein) for a spatially homogeneous mag- 
netic field. 

It is widely accepted (see, e.g., Takatsuka & Tama- 
gaki 1993, and references therein), that neutrons and/or 
protons in the NS core can be superfluid. Let Tcp and 
Ten be the critical temperatures of the proton and neu- 
tron superfluidities, respectively. Theoretical values of Tcp 
and Ten are strongly model dependent and range from 
10^ to 10^" K. For T < Tc, the superfluidity of nucleons 
strongly suppresses the neutrino luminosity produced by 
the Urea and bremsstrahlung processes in the NS cores 
(e.g., Pethick 1992, Yakovlev & Levenfish 1995). This en- 
hances the relative importance of the neutrino synchrotron 
process. 

It is commonly thought that superfluid protons form 
the type II superconductor (see Sauls 1989, and references 
therein). Let us assume that the core of a newly born, hot 
and initially nonsuperfluid NS contains a magnetic field 
(e.g., primordial one, or generated at the NS formation 
stage). The transition to a superconducting state in the 
course of the stellar cooling would then be accompanied by 
a dramatic change in the spatial structure of the magnetic 
field. Initially homogeneous field splits into an ensemble of 
fluxoids which contain a superstrong magnetic field, em- 
bedded in the field-free superconducting medium. So far, 
the neutrino synchrotron radiation has been studied only 
for a homogeneous magnetic field. In the present paper, 
we show that after the superfiuidity onset, the neutrino 
synchrotron mechanism transforms into the neutrino pair 
emission due to the electron-fluxoid scattering. 

Physical conditions in the superconducting NS cores 
are described in Sect. 2. In Sect. 3, we present general for- 
malism of the electron-fluxoid (e/) scattering. The calcu- 
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lations of the neutrino emissivity are performed in Sect. 4. 

In Sect. 5, wc estimate the emissivities of two additional 
neutrino production mechanisms which are usually ne- 
glected: the bremsstrahlung processes due to electron- 
proton (ep) and electron-electron (ee) scattering. No cal- 
culation of these emissivities has been done earlier, to our 
knowledge, although the importance of the ep scattering 
has been mentioned by Schaab et al. (1996). In Sect. 6 
we show that the new mechanisms (especially ef and ee 
scattering) can be important in the NS cores with highly 
superfluid protons and neutrons since the superfluidity re- 
duces strongly the familiar neutrino generation reactions 
(Urea processes, nucleon-nucleon bremsstrahlung). 

2. Superconducting neutron stcir cores 

The density within a NS core is expected to range from 
about 0.5po at its outer edge (Lorenz et al. 1993), to (5- 
10)po at the star center. For densities p < 2po, matter con- 
sists of neutrons, with a few percent admixture of protons, 
electrons, and possibly muons. The fraction of muons is 
usually much smaller than that of electrons, and a simple 
npe model is a good approximation. At higher densities, 
P ^ 2/9o, other particles may appear, such as hyperons, 
condensed pions and/or kaons, or even free quarks. De- 
pendence of the neutrino emissivity on the composition of 
dense matter is reviewed by Pethick (1992). For simplicity, 
we will neglect exotic matter constituents, restricting our- 
selves to the simplest npe model. All constituents of the 
npe matter are degenerate. The electrons form an ultra- 
relativistic and almost ideal gas, with the electron Fermi 
energy /ig ?ic(37r^ne)^/^ ^ (100 300) MeV, where n,, is 
the electron number density. The neutrons and protons 
form strongly interacting Fermi liquids. The Fermi energy 
of neutrons is of the order of that of the electrons, while 
the Fermi energy of protons is much smaller. 

Both nucleon species, n and p, can be in a superfluid 
state (for a review, sec Takatsuka & Tamagaki 1993). The 
superfluidity is widely accepted to be of the BCS-type; the 
nn and/or pairing occurs due to nuclear attraction. Let 
us focus on the proton superfluidity. In view of a relatively 
low proton number density, the inter-proton distance is 
rather high, and the pp interaction is attractive in the ^Sq 
state. Thus, the proton pairing is most likely to be in the 

state. Recent reviews on the properties of supercon- 
ducting protons in the NS cores are given by Sauls (1989) 
and Bhattacharya & Srinivasan (1995). The proton super- 
conductivity is characterized by an isotropic energy gap 
in the single- (quasi) particle spectrum Ap, which depends 
on temperature and density. The critical temperature Tcp 
is related to the zero-temperature energy gap Aj,(0) by 

= Ap(0)/(1.76fcB) (Lifshitz & Pitaevskii 1980), where 
fee is the Boltzmann constant. An important parameter of 
superconducting protons is the pp correlation length ^ ; it 
measures the size of a pp Cooper pair. In the BCS model, 
^ is related to Ap and the proton Fermi velocity vpp by 



^ = hvFp/{'n-Ap). The zero temperature value of ^ will be 
denoted by • Another important parameter is the pene- 
tration depth A of the magnetic field into a proton super- 
conductor. In the case of A ^ (which corresponds to the 
conditions prevailing in the NS cores), the zero tempera- 
ture penetration depth, Aq, is determined by the proton 
plasma frequency, cOp (e.g., Lifshitz & Pitaevskii 1980), 



Ao = 



1/2 



(1) 



where rip is the proton number density, and m* is the 
proton effective mass, that can differ from the bare mass 
nip due to the many-body effects in dense matter. Simple 
estimates yield typical values of of a few fm, and Aq of 
a few ten fm. If so. Aq ^ ^o. which means that protons 
constitute a type II superconductor. Notice that the values 
of Ao and are model dependent, and one cannot exclude 
the case of Aq ^ Co . in which the proton superconductivity 
could be of type I, although we will not consider this case 
in the present article (some comments on the case of a 
type I proton superconductor are given at the end of Sect. 
6). 

For T Tcp, both ^ and A diverge as (Tcp - T)-V2, 
while for T <^ Tcp they can be replaced by their zero tem- 
perature values, Co and Aq. However, we need to know the 
temperature dependence of A for all temperatures below 
Tcp. In what follows, we will approximate its tempera- 
ture dependence by the Gorter-Casimir formula (Tilley & 
Tilley 1990), 



A: 



An 



VI - {T/TcpY 



(2) 



The transition to the type II superconductivity dur- 
ing the NS cooling is accompanied by the formation of 
quantized flux tubes (Abrikosov fluxoids), parallel to the 
initial local magnetic field B. Each fluxoid carries an ele- 
mentary magnetic flux 0o = nhc/e. The number of flux- 
oids per unit area perpendicular to the initial field is 
A/f = B/(l)o, and the mean distance between the flux- 
oids is dp = [2(j)o/{\/3B)]^^^. Simple estimate yields dp ~ 
1500/(i3i3)^/2 fm, where B13 = B/{10^^ G). A fluxoid has 
a small central core of radius ~ C containing normal pro- 
tons. A typical fluxoid radius is A. Just after the supercon- 
ductivity onset (T < Tcp), this radius is large and the flux- 
oids fill all the available space. When temperature drops 
to about O.STcp, A reduces nearly to its zero-temperature 
value Aq. Thus, A becomes much smaller than the inter- 
fluxoid distance dp, for the magnetic fields B < lO"'^"'' G 
which we will consider in the present article. The maxi- 
mum value of B is reached at the fluxoid axis, Bmax — 
[(/)o/(7rA^)] log(A/C). In our case A >• C, and the magnetic 
field profile atr^^is given by (e.g., Lifshitz & Pitaevskii 
1980): 



B{r) 



27rA2 



(3) 
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where Ko{x) is a McDonald function. In particular, for 
r > A, one has B{r) « (/)o(87rrA^)"^/^ exp(-r/A). 

The superconducting state is destroyed when dp < i^, 
which corresponds to magnetic fields i? > Bc2 — (f>o / i'^^'^) , 
whose typical values are ^ 10^^ G. We do not consider 
such strong fields. 

Let us mention that the fluxoids may migrate slowly 
outward the stellar core due to the buoyancy forces (Mus- 
limov & Tsygan 1985, Jones 1987, Srinivasan et al. 1990). 

3. General formalism 

Consider the neutrino-pair emission due to scattering 
of strongly degenerate, relativistic electrons off fluxoids 
(Sect. 2), 



f 



(4) 



We will treat this process using the standard perturba- 
tion theory with free electrons in nonperturbed states. 
The process (^ is similar to the well known neutrino-pair 
bremsstrahlung due to scattering of electrons by atomic 
nuclei (see, e.g., Festa and Ruderman 1969, Soyeur and 
Brown 1979, Haensel et al. 1996, and references therein). 
It is described by two second-order diagrams, where one 
(electromagnetic) vertex is associated with electron scat- 
tering by the fluxoid magnetic fields, while the other (four- 
tail) vertex is due to the neutrino-pair emission. 

In what follows, we will mainly use the units in which 
c = ?i = fee = 1 although we will insert the ordinary 
physical units in the final expressions. 

The neutrino energy loss rate (emissivity) Qaux (ergs 
cm a ^^) of process (H) can be written as 



Qfiux — 



(2^) 



10 



dp J dp' J dKj dK 



X S{e~e' ^Lu)S{p,~p',-k,)L,f{l- f)W, (5) 

where A/p is the fluxoid surface number density defined 
in Sect. 2, P = (e,p) and P' = {s',p') are, respec- 
tively, the electron 4-momenta in the initial and final 
states, K„ — {LJi,,'ki,) and — (ujl^'kl) are the 4- 
momenta of the neutrino and of the antineutrino, and 
K = K„ + Kl = (w, k) is the 4- momentum of the neutrino 
pair (lu = Lu,^ + lu'^ and k = + k'^^). Furthermore, 



/ 



1 + cxp 



£ — /i 

T 



(6) 



is the Fermi-Dirac function for the initial electron state, 
and /' = f{e') is the same function for the final electron 
state. The ^-functions describe energy conservation and 
momentum conservation along the fiuxoid axis (the axis 
z). Finally, W is the differential transition rate 



(^2 



1 



2 (2w,)(20(2e)(2e') 



E 



Ml- 



(7) 



where Gf = 1.436 x 10^"'^ erg cm"^ is the Fermi weak 
interaction constant, and |Mp is the squared transition 
matrix element. Summation is over the electron spin states 
a before and after scattering and over the neutrino flavors 
(i/e, Pfi,, Vt)- The neutrino energies are assumed to be much 
lower than the intermediate boson mass (^--^80 GeV) . Then 
the standard approach yields 

J2 = y dq5(p - q - p' - k)AM^* 



X Tr(i^^O"i^^O^) 
X Tr Lf}j{P' + me)La,i{P + Tne) 

= r„G(F - Q)7,; -I- 7,G(P' + Q)r„ 
P + m. 



G{P) 



p2 



a _,5 



(8) 

(9) 
(10) 



Here, q = p — p' — k is a momentum transfer (in the {x, y)- 
plane) from an electron to a fluxoid, Q^P — P' — K = 
(f2, q) is an appropriate 4-momentum transfer (with no 
energy transfer, U — 0), A = A(q) is a 2-dimensional 
Fourier transform of the magnetic-field vector-potential, 
which lies in the {x, y) plane, and me is the electron mass. 
Greek indices a and P run over (0,1,2,3) and Latin ones i 
and j refer to the spatial components (1,2,3); G(P) is the 
free-electron propagator, 7" is a Dirac matrix, upper bar 
denotes Dirac conjugate, and P = Pq,7" (Berestetskii et 
al. 1982). Furthermore, Gy and Ca are the vector and the 
axial vector weak interaction constants, respectively. For 
the emission of the electron neutrinos (charged + neutral 
currents), one has Gy = 2sin^0w + 0.5 and Ca — 0.5, 
while for the emission of the muonic or the tauonic neu- 
trinos (neutral currents only), Cy = 2 sin^ 6'w — 0.5 and 
'w is the Weinberg angle, sin^ 6'w — 0.23. 
Using the identity (Berestetskii et al. 1982) 



C'^ = -0.5 



/ dk. / dk; 5(4) {K-K,- K) ^ 
J J ^i^^l 
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we obtain 



Qfiv 



''ClAfF 



12(27r)" 



dp Jdp' JdkSip, 



ee 



X Tr [(/" + me)i„^(F + me)Z^J 



(12) 



(13) 



where g""^ = diag(l, —1, —1, —1) is the metric tensor. The 
integration in Eq. (12) is to be carried out over the domain 
in which > 0. Notice that Eq. (12) reproduces the well- 
known expression for the neutrino-pair bremsstrahlung of 
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electrons which scatter off atomic nuclei (e.g., Haensel 
et al. 1996). For this purpose, one needs only to replace 
27rA/Fe^(5(pz - p'^ - kz)A'A^* ni\U{q)\'^Sij5io, where ni 
is the number density of nuclei and J7(q) is the Fourier 
transform of the electron-nucleus Coulomb potential. 

The 2-dimensional Fourier transform B{q) = Bz{q) of 
the fl-uxoid magnetic field (ph in cylindrical coordinates is 



(14) 



where Jo (a;) is a Bessel function and go = l/-^- Then, 
using cylindrical gauge, we have A(q) = ~iej^ B{q)/q, 

= (B X q)/{Bq) being a unit vector. Accordingly, in 
Eq. (12) A'A^* = {B{q)lqf e\e\, where i, j = 1 or 2, for 
nonvanishing components. 

Equations (12) and (13) determine the neutrino emis- 
sivity for any degree of electron degeneracy and relativism. 
We are interested in the case of ultrarelativistic, strongly 
degenerate electrons (Sect. 2). In the relativistic limit, te- 
dious but straightforward calculations yield: 



(15) 



J4 



g^[2(eAp)(e^p')+ee'-pp'] + Sif^ 



uv 



(ee' - pp') [K^ - 2(e^k)2] 



2if4[2(e^p)(e^p')+ee'-pp'] {- -- 

\u V 



8K' 



uv 



(qk)(e^k)2 



1 1 



(2qk - K^) 



A T^ii 1 ^ ( gap' 
V u t; 



(16) 



where we introduce the notations 
u = ]^{P' + Kf = P' K +]^K\ 

1 , x2 1 9 

V = --[P-Kf =PK--K\ 
Using Eqs. (12), (|l5|), and (16), we get 

^AAf / dp y dp'! A\^5{p,-p',-k,) 



(17) 



12(27r)E 



(18) 



For practical applications, it is convenient to express 
Qflux in the form 



2.66 X 10^^ Bi3 TI 



Ip I lip 



HQ m 



1/2 



T 



1/2 



L erg s ^ cm 



(19) 



Here, no = 0.16 fm^^ is the standard (nuclear) number 
density, Tg is temperature in the units of 10^ K, and Bi^ 
is defined in Sect. 2. Numerical expression for Qaux takes 
into account emission of J/g, Vp, v^, so that C\ w 1.675. 
In Eq. (19), we introduced the dimensionless quantity L, 
defined by 



L 



189 



(27r)9T6g, 



-j dpdp'dk(5(p^ -p'^ - fc^ 



/ ££' 



(20) 



with J+ given by Eq. (16). Notice that Eq. (20) is vahd 
for any axially symmetric distribution of the fluxoid mag- 
netic field B{r) although we will use a specific distribution, 
Eqs. (^ and (14). An analogous quantity L was introduced 
in the case of the neutrino-pair bremsstrahlung due to the 
electron-nucleus scattering. In that case, it had a meaning 
of the Coulomb logarithm (Haensel et al. 1996). 

4. Practical evaluation of L 

Let us calculate L from Eq. (20). Since the electrons are 
strongly degenerate (Sect. 2), the main contribution to L 
comes from those electron transitions, in which the elec- 
tron momenta p and p' lie in the narrow thermal shell 
around the Fermi surface, |e — /i| < T and je' — /x| < T. 
In Eq. (20) we may set dp = e'^ de dfi, where dil is a solid 
angle element in the direction of p. We also put e = p 
since the electrons are ultrarelativistic (Sect. 2). Then we 
can set e = p^e in all smooth functions of s. Afterwards, 
the integration over e is standard: 



J defil-n 



(21) 



The integrand in Eq. (20) depends obviously on the rela- 
tive azimuthal positions of p, p' and k. Thus we can place 
p in the (a:z)-plane without any loss of generality. 

We shall mainly consider the case of qo <^ ppc typi- 
cal for the NS cores (Sect. 2). It is convenient to replace 
the integration over p' in Eq. (20) by the integration over 
q. The integration is simplified because the main contri- 
bution comes from the values of g <C ppe- Thus we can 
use the approximation of small-angle scattering. Since k 
can be comparable to q, we should keep two first terms 
in Eq. (16). Using the inequalities k <^ p and q ^ p we 
obtain 



J4 



8(a;2 



(22) 
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where we have neglected the quantity ee' — pp' = [(q + 4-1- Low-temperature, bremsstrahlung regime, T <C Tq 



^]/2 in the first term of (16), and have used the 



approximate expressions u !v v ^ pgr, which foUow from 
Eqs. (17). The subscript r denotes the vector component 
along P', Qr = Qx sin 9, where 6 is an angle between p and 
the z-axis (electron pitch- angle). 

Using energy — momentum conservation, in analogy 
with the results of Haensel et al. (1996) we obtain 



2 7 2 

qr{2uj - 



(23) 



where kt is a component of k transverse to p. 

The condition lo^ > k^ requires > 0, w > <7r/2 and 
kf < kQ — qr{2uj — qr)- Taking into account Eq. (|2l|), we 
obtain 



L = 



189 



4>Q qr 



9,72 



_ 1 



Akt h {kl-ki). 



(24) 



The integrals converge rapidly due to sharp decrease of 
B[q), Eq. (14). In view of this, we could extend the inte- 
gration over qx and qy up to cxd, and the lower integration 
limit over qy down to —00. 
Integrating over kt , we have 

189 , , 

dqx / aqy 



L 



B{q) 



de SIT? 

Jo 



duj 

90 J 

which can be transformed as 
189 g3_ r 
' lo 



(25) 



L = 



327r6 re 
duj 



de sir? 



00 



9r./2 



{qi 

2 



X 



qlY" 



(26) 



These equations take into account weak inelastic- 
ity of the e/ scattering, in analogy with the neutrino 
bremsstrahlung due to the electron — nucleus collisions 
(Haensel et al. 1996). 

As can be seen from Eq. (26), L depends on the dimen- 
sionless parameter t^ ~ T/Tq: 

* \ 1/2 r / _ N 41 -1/2 

TO* 710 ^ ' 



0.00786 Tc 



1 - 



T 



T 



(27) 



where Tq = Tp^l - [T/T^pY, and Tp = hujp/kB is the 
proton plasma temperature corresponding to the proton 
plasma frequency ujp (Sect. 2). Typically, Tp ^ I MeV, 
under the conditions in the NS cores. Let us analyze two 
extreme cases: T < To and T > To. 



In the limiting case of to ^ 1; Eq- (26) gives L — 7r/4. 
Then the neutrino emissivity is then given by 



Qflux — 



9072 f?c^ 



{k^TfqaUF 



2.09 X lO^^BiaTe 



no TO* 



1/2 



1 - 



T 

t7 



cp 



1/2 



erg s ^ cm ^. 



(28) 



The neutrino emission in the low-temperature regime 
is very similar to that due to the electron-nucleus 
bremsstrahlung (e.g., Haensel et al. 1996). In particular, 
the emissivity Qaux is proportional to T®. Therefore we 
will call this regime the bremsstrahlung regime. 

For the sake of completeness, we have also examined a 
not too realistic case of qa ^ ppe at T <C To . The analysis 
is based on Eq. (20). In this case, the integration can be 
simplified, and reduces to 



2/0 



L - I d^j 



1 



2y' 



1-2/' 



■ log(2/) 



(29) 



where y = q/{2p-pc)- Now, T is a function of the only 
dimensionless parameter (see Sect. 2) 



yo 



2pFe \37r?ic rripC 



1/2 



T 

TcT) 



1/2 / X 1/6 

0.0165 ( ^ 



p 



no 



T 

t7 



cp 



1/2 



-I 1/2 



(30) 



In the limit of ?/o ^ 1 we reproduce our basic result L = 
7r/4. For a finite yo, the function L given by Eq. ( p9| ) can be 
fitted (with relative error < 2.4 %) by a simple expression 



L 



4^/1 + 0.70572/2- 



(31) 



4-.2. High-temperature, synchrotron regime, T ^ To 



In the case of to ^ Ij from Eq. (26) we obtain an asymp- 
totic form L « (189/7re) C(5)/to, where C(5) « 1.037 is the 
value of the Riemann zeta-function. This yields 



C(5)G|e202c2 



127r6 r c7 



■ 90 (fceT) Wf « 6.89 x 10^^ B 



13 



X TS 



no TO* 



T 

J- r.r] 



erg s ^ cm ^. 



(32) 



The temperature dependence (Q cx T^) is the same as for 
the synchrotron emission of neutrinos by electrons (e.g.. 
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Kaminker et al. 1991). Thus we will refer to the high- 
temperature regime as the synchrotron regime. 

Let us compare the e/-scattering emissivity (32) with 
the emissivity Qsyn of the "purely synchrotron process" 
for the most important case of ^ T ^ M^^s ( Me is 
the electron chemical potential and LUg = eBc/ fj,e is the 
electron gyrofrequency) . For the uniform magnetic field 
Kaminker et al. (1991) obtained: 



isyn 



2C(5) e^GlClikBT)'' 
97r5 



9.04 X lO^^Bf^T^ 



ergs cm 



(33) 



(The numerical factor 9.04 differs from a factor 8.97 in 
Eq. (17) of Kaminker et al. 1991, because in the present 
paper we use more accurate value of the Fermi constant 
Gf.) 

We start with a brief discussion of the e/ scattering at 
to ^ 1 for an arbitrary distribution of the fluxoid mag- 
netic field B{q). In this case, we can set ^ in Eq. (25). 
Then the integrations over 9 and lo are decoupled from 
the integration over q, and can be done analytically. This 
yields the emissivity 



Qtiu.= j^C{5)GlClT'UF I dqi?2(q), 



(34) 



where the integration should be done over the entire 
(feiQy) plane. The emissivity is seen to be indepen- 
dent of the electron Fermi momentum. In particular, in 
the limit of a (quasi) uniform magnetic field B{r) — 
Bo, B{q) = 47r2 Bo (5*^^(q), we can replace B{q)'^ 
47r2 Bl (5(2) (q) /AAp , and obtain 



C(5) e^GlClikBTf Bl 
37r5 



Comparing ( p3^ ) and (33), we have Qsyn = (2/3) Qaux- 
The same factor 2/3 is obtained, if we treat B in Qsyn as 
a local magnetic field of a fluxoid, average Qsyn oc B'^ 
over a lattice of fluxoids with the magnetic field (^), and 
compare the result with Qaux, Eq. (32). Let us empha- 
size, that this emissivity Qsyn, averaged over a nonuniform 
magnetic field, can differ from the synchrotron emissivity 
Qsyn in the initial uniform field -B by a large factor of 
52/52 _ 0gg2/(47riJ) ^ (dF/A)2, where dp is an inter- 
fluxoid distance (Sect. 2). This increase of Qsyn is caused 
by the magnetic field enhancement within the fluxoids due 
to magnetic flux conservation. The enhancement is much 
stronger than the reduction of neutrino-emission space oc- 
cupied by the fluxoids. 

A detailed analysis shows, that the difference by a 
factor of 2/3 comes from momentum space available for 
neutrino-pair momentum k. The space is different in the 
case of synchrotron radiation in a uniform magnetic field 
and in the case of ef scattering by magnetic inhomo- 
geneities (in our case — fluxoids). Let us fix the electron 



pitch angle and average over the electron Larmor rota- 
tion. In the both cases k is concentrated to the cone of 
a given 9 (neutrinos are emitted predominantly within a 
narrow angle interval 59 in the direction of the electron 
pitch angle 9), and typically, fc ~ T. In the case of the 
ef scattering, our estimates yield 50flux ~ a/ Qo/T. The 
synchrotron radiation consists of many discrete cyclotron 
harmonics s with typical values s ~ T/cog ^ 1, for the 
quasiclassical regime of interest. For a fixed 9, the radia- 
tion in each harmonics is peaked at the pitch-angle cone, 
and maximum neutrino momenta k are restricted by a 
parabolic surface (Kaminker et al. 1991). The maximum 
value of k in the pitch-angle direction is slu g / {sin 9)'^ , and, 
typically, 50syn ~ s~'^^^. Replacing a sum over harmonics 
by an integral we obtain a smoothed synchrotron k space. 
It fills only 2/3 of the available momentum space, due to 
the parabolic momentum space restriction for each har- 
monics (e.g., the surface area below the parabolic curve 
y = l — x2at0<a;<lis exactly 2/3 of the surface 
area below a straight segment y = 1). A transition from 
the synchrotron formula (33) to the fluxoid formula (32) 
is expected to occur at (56'flux ^ 59syn- This is equivalent 
to the condition A ~ {c/u;g)s~^^^ , which takes place when 
the scale of the magnetic field variation A becomes compa- 
rable to a distance along which a neutrino pair is emitted 
in a uniform magnetic field, i.e., to the s^^^-th part of the 
Larmor radius c/ujg. 

4-3. Overall fit and overview 

Combining the results of Sects. 4.1 and 4.2, we can pro- 
pose the following fitting expression for the quantity L in 
Eq. (19) at yo < 1: 



(35) L^LoUV, 



(36) 



in = - 



TT 0.260to-f 0.0133 



4 t^-f 0.25 to + 0.0133' 



U = 



27 + 1 
37 + 1' 



.38 X 10~i2 TgBf^ 



V = 1 



^y3/6 .^,3/2 

ne J \mp 
AttBX^ 



T 

t7 



cp 



-3/2 



(37) 



00 



m„ no 
1 + 0.00210 Bi3 ^ 



T 



n -1 



(38) 



Here, Lq is the analytic expression which fits the results 
of our numerical calculations of L from Eq. (26) with er- 
ror < 1 % for any to- The factor U ensures, somewhat 
arbitrarily, the difference by a factor of 2/3 between the 
cases of weakly and strongly nonuniform magnetic fields 
(Sect. 4.2). The factor V provides a smooth transition 
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from the ef scattering to the pure synchrotron emission 
in the case when a fiuxoid radius A becomes very large 
(T Tcp), the neighboring fluxoids overlap and the over- 
all magnetic field becomes nearly uniform. 

Using our results, we can follow the evolution of Qaux 
in the course of the superconductivity onset. If T > Tcp, 
the emissivity Qflux is given by Eq. (33) since it is es- 
sentially the same as the synchrotron emissivity Q^yn = 
QsyniB) in a locally uniform 'primordial' magnetic field 
B = B. After T falls only slightly below Tcp, the factor U 
transforms from U ~ 2/3 ('pure synchrotron') to J7 = 1 
(e/ scattering in the synchrotron regime). The fluxoid 
structure is still not very pronounced, i.e., V w AttBX^ / (f>Q, 
Qflux ~ Qsyn. When T decreases to about 0.8 T^p, the flux- 
oid radius A nearly achieves its zero-temperature value Aq, 
and V ~ I. Accordingly, Qflux/Qsyn grows by a factor of 
(o?f/Ao)^ due to the magnetic field confinement within the 
fluxoids (Sect. 4.2). Simultaneously, <o decreases from very 
large values at T ^ T^p to to — T^p/Tp. 

Let us first consider the case in which the proton 
plasma temperature is Tp <^ Tcp. Then the emissiv- 
ity Qflux is enhanced with respect to Qsyn by a factor 
of (dp/Ao)^ at T « O.STcp, and this enhancement re- 
mains nearly constant over a wide temperature range 
down to Tp. Within this range, the ef scattering oper- 
ates in the synchrotron regime. At lower temperatures, 
T < Tp, the synchrotron regime transforms into the 
bremsstrahlung regime (Sect. 4.1) and we have Qaux ~ 
Qiyiidp/Xo)^ (T/Tp). Thus, for T < Tp the emissivity 
Qflux decreases with respect to Qsyn, and may become 
lower than Q^yn. 

In the opposite case of Tp > Tcp, the ef scatter- 
ing operates in the synchrotron regime only in a narrow 
temperature range below Tcp, and transforms into the 
bremsstrahlung regime at lower T. In the latter regime, 
we have Qhux ~ Qi%idF/X)^iT/Tp). 

5. Neutrino bremsstrahlung due to ep and ee scat- 
tering 

As will be shown in Sect. 6, the ef scattering mechanism 
can be important, if protons and neutrons in a NS core 
are highly superfluid (have high critical temperatures Tcp, 
Ten), SO that the rates of the "standard" Urea and nucleon 
bremsstrahlung processes are strongly suppressed. Under 
such conditions, one should carefully take into account 
all the neutrino processes, even if they are negligible for 
relatively low critical temperatures. 

In this section, we consider two additional neutrino 
production mechanisms, the neutrino-pair bremsstrahlung 
due to the ep and ee scattering: 

e+p-^e + p + iy + D, 

e + e~*e + e + u + D. (39) 



As far as we know, these mechanisms were neglected in 
the neutron-star cooling simulations. 

Let us restrict ourselves to simple estimates of the 
emissivities Qep and Qee', detailed analysis will be pub- 
lished elsewhere (Kaminker et al. in preparation). We 
start from the well-known emissivity Qez of the electron- 
nucleus bremsstrahlung (e.g., Haensel et al. 1996). In the 
particular case of the ep collisions of relativistic, strongly 
degenerate electrons in a nonsuperfluid matter, we obtain 

Qi'J- 'l^Jllf {kBTfnpLcp, (40) 

where Lep is the Coulomb logarithm. The above formula 
has been derived assuming, that an energy transfer to a 
proton in a neutrino-emission act is much smaller than 
T. As applied to our case of strongly degenerate protons, 
the latter means that the Pauli principle imposes no re- 
strictions on the states of protons. This is true as long as 
T ^ (7s, where Qs is an inverse length of plasma screening 
of the Coulomb ep interaction (see below). However, in 
reality we have just an opposite regime, T <C gs, and Eq. 
(|0|) does not hold. 

To estimate Qep at T <C (/s, we assume that Qep oc Vep, 
where Vep is some effective ep collision frequency. If our as- 
sumption is correct, we have Qep ~ Qlp'' i^ep/i^ep , where 
i/ip^ and Vep refer to the high- and low- temperature 
cases, respectively. For high temperatures, we will use the 
effective frequency of nearly elastic collisions vip'' (e.g., 
Yakovlev and Urpin 1980) which determines the electric 
and thermal conductivities of degenerate electrons. For 
T <^ Qs, we employ the collision frequency Vep that de- 
termines the electron thermal conductivity (e.g., Gnedin 
and Yakovlev 1995). Performing the above rescaling, we 
obtain 

_ 3.42 X 10^^ fm^V f^V'^ 

Vl \mpj \npj 
X RepT^ ergscm-3s"\ (41) 

where j/s = ?i9s/(2pFc) is the plasma screening parameter. 
We additionally introduce the factor Rep that describes 
suppression of the emissivity by the proton superfluid- 
ity. According to Gnedin & Yakovlev (1995) the plasma 
screening in npe matter of the NS cores is defined by 

y!-^(^ + ^Zp), (42) 

where m* — /ie/c^. The first term in brackets comes from 
the electron screening while the second is due to the pro- 
ton screening. The latter contains the factor Zp which de- 
scribes reduction of the proton screening by the proton 
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superfluidity. According to Gnedin & Yakovlev (1995) this 
factor is fitted as 



1/2 



= (^0.9443 + v/(0.0557)2 + (0.1886i; )2j 
X exp (1.753 - \/(l-753)2 + i;2^ , 
wfiere v is tfie proton superfluid gap parameter 



(43) 



= Wl- 




kBT 



X I 1.456 - 0.157 \/ ^ + 1.764 ^ 



If protons are normal (T > T^p), one has R, 



(44) 



"-ep 



in Eqs. (41) and ([12[). The factor Rep, which damps Qep 
at T < Tcp, should be the same as the factor i?„p that 
describes reduction of the neutrino emission in np colli- 
sions by the proton superfluidity (Yakovlev and Levenfish 
1995): 



Re 



1 



2.732 



a exp (1.3O6 - V(l-306)2 +w2^ 
1.732 6^ exp (^3.303- \/(3.303)2 + 4w 
v/(0.0018)2 + (0.3815w)2, 



(45) 



with a = 0.9982 
5 = 0.3949 + v/(0.6051)2 + (0.2666w)2. 

Now consider the neutrino emissivity Qee due to ee 
scattering. Under the same assumption, we obtain Qee — 
Qepi^ee/i^ep, wlicrc Vee IS the effective ee collision fre- 
quency that determines the electron thermal conductivity 
(e.g., Gnedin & Yakovlev 1995). Then 



Qe 



378?i^ciOy3p|, 



1.07 X 10^4 



He {kBTf 
X 1/3 

no 



ergs cm ^ s ^. 



(46) 



Although we do not calculate Qpe and Qee exactly, we 
believe that Eqs. (41) and (46) give reliable estimates. The 
emissivity Qep is reduced exponentially by the strong pro- 
ton superfluidity, while Qee is affected by the superfluidity 
in a much weaker manner, only through the plasma screen- 
ing parameter ( ^2|) . If protons are normal they provide the 
major contribution into the plasma screening. If they are 
strongly superfluid, a weaker electron screening becomes 
important, which enhances Qee, but not to a great extent 
(see below). 

6. Discussion 

Figures 1 and 2 illustrate the efficiency of various neutrino 
production mechanisms in npe matter of the NS cores. We 
compare the neutrino emissivity due to the e/ scatter- 
ing or due to the electron synchrotron radiation (Sect. 4) 



with those produced by the neutron and proton branches 
of the modified Urea reactions (Friman & Maxwell 1979, 
Yakovlev & Levenfish 1995), 

71 + 71— s-n + p + e + Pe, n + p + e n + n + h'e ■ 
n + p ^ p + p + e + De, p + p + e ^ n + p + i^e ; (47) 

by the nucleon-nuclcon bremsstrahlung processes (Friman 
& MaxweU 1979, Yakovlev & Levenfish 1995) 

n + n^n + n + v + D, 
n + p^ji + p + v + D, 

p + p^p + p + v + D] (48) 

and also by the ep and ee bremsstrahlung processes (39) 
considered in Sect. 5. For better presentation, the total 
emissivity of the both branches of the modified Urea re- 
actions (47) is represented by a single curve 'Murca\ 
and the total emissivity from the three nucleon-nucleon 
bremsstrahlung reactions (48) are represented by a sin- 
gle curve 'NN'. We adopted a moderately stiff equation of 
state of Prakash et al. (1988) (the same version as used 
by Page & Applegate 1992) and choose matter density 
p = 5.6 X IQ^^ g cm-3 (7ip = ne~ 0.0207 fm'^, n„ « 0.323 
fm"'^), well below the threshold density 1.30 x lO""^^ g cm~'^ 
at which the direct Urea process becomes allowed (Lat- 
timer et al. 1991). Therefore our curves correspond to the 
standard neutrino emission (no exotic cooling agents, see, 
e.g., Pethick 1992). We show temperature dependence of 
various neutrino emissivities in the range from 10^ K to 
5 X 10^ K, which is the most important for the cooling 
theories of NS. The nucleon effective masses are set equal 
to 0.7 of the masses of bare particles. The neutrino emis- 
sivities vary slowly with density (below the direct Urea 
threshold), i.e., we present a typical situation within the 
NS cores. In both figures, we assume the presence of the 
'primordial' magnetic field B = 10^"^, 10^^ or lO^'' G. 

Figure 1 corresponds to a nonsuperfluid matter. If the 
magnetic field is absent, the Murca processes are seen to 
be the most important ones, in the displayed temperature 
range. The NN bremsstrahlung is much weaker, and the 
ep and ee bremsstrahlung processes are even weaker. The 
emissivities of all these mechanisms vary as T^, whereas 
the emissivity of the synchrotron radiation varies as T^. 
Thus the synchrotron radiation is negligible for T > 10^ 
K but becomes more important with decreasing T. At 
T ~ 10* K and B = 10^^ G the synchrotron dominates 
over all other mechanisms except the Murca processes, 
and at slightly lower T it becomes the most efficient of 
all mechanisms. However, such temperatures are too low 
to be important for the cooling theories: the neutrino lu- 
minosity of NS becomes then smaller than the photon lu- 
minosity from the stellar surface, and the neutrino emis- 
sion looses its significance as a source of the NS cooling. 
Figure 2 shows how the neutrino emission displayed in 
Figure 1 is modified by the nucleon superfluidity. In our 
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Fig. 1. Neutrino energy loss rates from various processes, ver- 
sus temperature in a non-superfluid npe matter of density 
5.6 X lO'^* g cm~^. Murca — total contribution of two branches 
of the modifiod Urea process, Eq. (47); NN total contribu- 
tion of three nucleon bremsstrahlung processes, Eq. (48); ep 
and ee — bremsstrahlung processes, Eq. (39), due to the ep 
and ee scattering, respectively; solid lines — synchrotron radi- 
ation of electrons for three values of the magnetic field. 



example, the protons and neutrons become superfiuid at 
Tcp = Ten = 2 X 10^ K. These are typical, moderate 
critical temperatures for the nucleon superfluidity (e.g., 
Takatsuka & Tamagaki 1993). The factors which describe 
superfiuid suppression of the modified Urea processes (47) 
and NN collisions (48) are taken from Yakovlev & Lev- 
enfish (1995). At T < 10^ K the nucleon superfluidity is 
seen to strongly reduce all the neutrino production pro- 
cesses which involve nucleons. If the magnetic field were 
absent, the ee scattering would be the main neutrino gen- 
eration mechanism at T < 5 x 10^ K. However, the ef 
scattering in the presence of 5 > 10""^^ G also becomes 
important, and it can even dominate. After the superfluid- 
ity onset (T < Tcp), the magnetic field splits into fiuxoids 
and the synchrotron radiation transforms into the radi- 



Fig. 2. Same as in Fig. 1 but in the presence of proton and 
neutron superfluidities with Tcp — Ten — 2 X 10" K. At T < 
the initially uniform magnetic field splits into fiuxoids, and the 

synchrotron emission transforms into the emission due to the 
ef scattering. The nucleon superfluidity suppresses all mecha- 
nisms but the ee and ef scattering. 

ation due to the ef scattering. This enhances Qqux ow- 
ing to the field amplification within the fiuxoids (Sect. 4). 
The enhancement is more pronounced at lower B ~ 10^^ 
G at which the field amplification is stronger. This is a 
rare case in which the neutrino emissivity increases with 
decreasing T. For the conditions displayed in Fig. 2, the 
proton plasma temperature Tp = 5.47 x 10^" K is much 
higher than the superfiuid critical temperature Tcp. There- 
fore, the synchrotron regime (Sect. 4.2) in the ef scatter- 
ing operates only at T < Tcp, during the phase of the 
fiuxoid formation. Very soon after T falls below Tcp, the 
synchrotron regime transforms into the bremsstrahlung 
regime (Sect. 4.1) which operates further with decreasing 
T. 

We have analysed a number of cases, varying the pro- 
ton and neutron critical temperatures Tcp and Tc„. Our 
principal conclusion is that the standard neutrino pro- 
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duction mechanisms (47) and (48) dominate in the tem- 
perature domain of practical interest if either Tcp and/or 
Ten are not too high (not higher than about 10^ K). If, 
however, both critical temperatures are high, the situa- 
tion is similar to that shown in Fig. 2. The superfluidity 
suppresses all the traditional neutrino generation mecha- 
nisms, and the main neutrino production at T < 5 x 10^ 
K comes either from the ef or the ee scattering. 

In principle, neutrinos can be generated also in the pp 
collisions of normal protons in the non-superfluid cores 
of the fluxoids (Sect. 2) as well as in the en collisions. 
However, the non-superfluid cores are very thin, the emis- 
sion volume is minor, and the first process is inefficient. 
The neutrino bremsstrahlung due to en scattering is also 
inefficient since it occurs through electromagnetic interac- 
tion involving neutron magnetic moment (e.g., Baym et 
al. 1969) and since it is suppressed by the neutron super- 
fluidity. 

While it is commonly accepted that protons form a 
type II superconductor, a possibility that they actually 
form a type I superconductor cannot be excluded. This 
uncertainty results from the lack of a precise knowledge 
of the nucleon-nucleon interaction in dense nuclear mat- 
ter, and from the approximations and deficiencies of the 
many-body theory of dense nucleon matter. In the case 
of a type I proton superconductor, which corresponds to 
^0 > "v/SAo (de Gennes 1966), cooling below Tcp is ex- 
pected to be accompanied by a transition of the magne- 
tized interior to an "intermediate state" (de Gennes 1966, 
Baym et al. 1969). The "intermediate state" would con- 
sist of alternating regions of normal matter containing 
magnetic flux, and superconducting regions exhibiting a 
complete expulsion of magnetic flux. The specific spatial 
structure of the "intermediate state" would result from 
the condition of the minimum of the thermodynamic po- 
tential at a fixed macroscopic magnetic fiux (de Gennes 
1966). Qualitatively, we expect that the presence of an 
"intermediate state" in the superconducting proton core 
will imply an enhancement of its electron synchrotron i?i> 
emissivity, compared to the case of a normal proton core, 
because of > {BY. However, for B < 10^"* G this en- 
hancement will be much smaller than that characteristic 
of the type II superconductor, in which magnetic field is 
confined to fluxoids. 

7. Conclusion 

We have considered (Sects. 3 and 4) the neutrino-pair ra- 
diation (||) due to scattering of relativistic, degenerate 
electrons off threads of quantum magnetic flux — fluxoids 
— in the npe matter within the superfluid NS cores. We 
have shown that this mechanism is a natural generaliza- 
tion of the synchrotron emission of neutrinos by electrons 
in a non-superconducting matter (with a locally uniform 
magnetic field) to the case in which the protons become su- 
perfluid and the magnetic field splits into the fluxoids. Ac- 



cording to our results, this mechanism can operate either 
in the synchrotron (Sect. 4.2) or in the bremsstrahlung 
(Sect. 4.1) regime. We have obtained a simple fitting ex- 
pression (^6|) which reproduces the emissivity (19) for any 
parameters of practical interest. 

Furthermore, we have estimated (Sect. 5) the neu- 
trino emissivities in two additional neutrino production 
mechanisms (39), the ep and ee bremsstrahlung processes. 
In a non-superfluid matter, these mechanisms are much 
weaker than the standard neutrino emission mechanisms, 
such as the modified Urea processes (47) and the nucleon- 
nucleon bremsstrahlung (48). If, however, both superfluid 
critical temperatures, Tcp and Ten, are higher than 10^ 
K, then the superfluidity strongly suppresses all the tra- 
ditional (standard) neutrino energy losses. In such a case, 
the main neutrino generation in a NS core at T < 5 x 10^ 
K occurs either via the electron-fluxoid or via the ee scat- 
tering (Sect. 6). Therefore our results can be of particular 
importance for simulating the cooling of NSs with highly 
superfluid cores. 
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